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Nematicity has emerged as a key feature of cuprate superconductors, but its link to other funda¬ 
mental properties such as superconductivity, charge order and the pseudogap remains unclear. Here 
we use measurements of transport anisotropy in YBa 2 Cu 30 a to distinguish two types of nematicity. 

The first is associated with short-range charge-density-wave modulations in a doping region near 
p = 0.12. It is detected in the Nernst coefficient, but not in the resistivity. The second type prevails 
at lower doping, where there are spin modulations but no charge modulations. In this case, the 
onset of in-plane anisotropy — detected in both the Nernst coefficient and the resistivity — follows 
a line in the temperature-doping phase diagram that tracks the pseudogap energy. We discuss two 
possible scenarios for the latter nematicity. 

PACS numbers: 74.72.Gh, 74.25.Dw, 74.25.F- 

I. INTRODUCTION 

It has become clear that charge-density-wave (GDW) 
order is a generic tendency of cuprate superconduc¬ 
tors [1-7]. In YBa 2 Cu 30 y (YBCO), short-range GDW 
modulations detected by X-ray diffraction (XRD) appear 
below a doping p ^ 0.16 [8,9]. Both the amplitude of the 
GDW modulations and their onset temperature, Txrd, 
peak at p = 0.12 (Fig. 1) [8,9]. NMR measurements 
in high magnetic fields [1,10] detect the abrupt onset of 
GDW order below a temperature Jnmr that also peaks 
at p = 0.12. This GDW order causes a reconstruction of 
the Fermi surface [11] detected as a sign change in the 
Hall [12,13] and Seebeck [14,15] coefficients, from positive 
at high temperature to negative at low temperature, due 
to the formation of a small electron-like pocket [16]. The 
temperature Tmax below which R}i{T) drops as a result of 
Fermi-surface reconstruction peaks at p = 0.12 [13]. The 
GDW order and associated reconstruction, both peaked 
at p = 0.12, cause a dip in Tc (ref. 17) (Fig. 1) and a lo¬ 
cal minimum in the superconducting upper critical field 
Hc 2 (ref. 18) at p = 0.12. 

The short-range GDW modulations are known to be 
anisotropic within the Cu02 planes of YBCO, [8,9] and 
they appear to break the rotational (C4) symmetry of 
the planes, [20,21] i.e. they have a nematic character. 

This should give rise to an in-plane anisotropy of trans¬ 
port. A large in-plane anisotropy of the Nernst coefficient 
v{T) has indeed been observed in YBCO (see Fig. 2) [22], 
and it does grow in tandem with the GDW modulations 
measured by XRD at p = 0.12 (see Fig. 3). 

In the original report [22], for dopings p = 0.12 and 
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higher (0.13,0.15 and 0.18) the Nernst anisotropy differ¬ 
ence was shown to collapse onto a common curve when 
plotted as {va — Vb) / T vs T / T* (Fig. 4 of ref. 22), 



FIG. 1: Temperature-doping phase diagram of YBCO, with 
the superconducting phase below Tc (black dots [17]) and 
the antiferromagnetic phase below Tn (green line). Short- 
range charge-density-wave (CDW) modulations are detected 
below Txrd (up triangles; [8] down triangles [9]). Tnem 
(red symbols) marks the onset of the rise in the resistivity 
anisotropy (circles; from ref. 19 and Fig. 5(b)) and in the 
Nernst anisotropy (squares; from Figs. 5(c) and 5(d)), as tem¬ 
perature is reduced. The blue and red lines are guides to the 
eye. They delineate the two regions of the phase diagram 
where two types of nematicity prevail — one associated with 
CDW modulations (blue), the other not (red). 
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FIG. 2: Anisotropy of the Nernst coefficient v{T) compared 
with the corresponding anisotropy of the conductivity (j{T) 
in YBCO, both plotted as ratios: vt, / Va (dots) and at, j Oa 
(curve) (adapted from ref. 22). The CuO chains that run 
along the b axis of the orthorhombic structure make a con¬ 
tribution to (Ti, such that ai, j (Ja > 1 at T = 300 K and 
CT(, / (Ja decreases below T = 150 K [19]. (a) At p = 0.078, 
both Vf, / i/a (dots) and ai, / Ua rise with decreasing tempera¬ 
ture, so that both reveal the emergence of nematicity in the 
Cu02 planes at low temperature, (b) At p = 0.12, nematic¬ 
ity is observed in ui, / Va, as a. huge anisotropy that rises with 
decreasing T. However, no corresponding anisotropy rise is 
detected in ai, j cra, which only shows a drop typical of that 
due to the CuO chains [19]. Those two comparisons therefore 
reveal an empirical difference in the nematicities at p = 0.078 
and at p = 0.12. 


where T* is the pseudogap temperature defined as the 
onset of the drop vuv^jT vs T below its constant value 
at high temperature (inset in Fig. 1 of ref. 22). (This def¬ 
inition was shown to be equivalent, within error bars, to 
the usual definition of T* from the downward deviation 
in the resistivity Pa{T) vs T from its linear dependence 
at high temperature (Fig. S4 of ref. 22).) What this col¬ 
lapse revealed is a common, very slow increase of Nernst 
anisotropy starting at T*, for p = 0.12 and higher (see 
Fig. 4(a)). This led to the conclusion that the pseudogap 
phase causes an extra anisotropy in transport. 

However, as shown in Fig. 4(b), the different curves 
scale just as well when plotted as T / Txrd, perhaps even 
better. Note, moreover, that the growth in (— i/f,) / T 
that occurs between T* and Txrd is only 1% of the 
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FIG. 3: Temperature dependence of the X-ray intensity as¬ 
sociated with the short-range CDW modulations in YBCO 
at a doping p = 0.12 and a magnetic held H — 15T (blue 
squares, left axis; from ref. 3). Note that above T = 66K 
the data at H — 0 are identical to the data at 7/ = 15T [3]. 
Txrd = 140 K marks the approximate onset of the gradual 
rise in CDW intensity (blue arrow). The in-plane anisotropy 
of the Nernst coefficient Vx{T) measured at the same doping 
and the same field is plotted as {va — Vb) ! T (red dots, right 
axis; from ref. 22). It is seen to track the CDW intensity very 
well, compelling evidence that the anisotropic CDW modu¬ 
lations cause the large Nernst anisotropy at that doping. At 
dopings above p = 0.12, the rise in the Nernst anisotropy 
continues to track the onset of CDW modulations, i.e. the 
different curves of {va — I'b) / T for different dopings scale 
when plotted as T /Txrd (see Fig. 4(b)). 


growth that takes place below Txrd , tracking the XRD 
intensity (Fig. 3). Therefore, the large Nernst anisotropy 
in YBCO at p = 0.12 (and higher doping) is more likely 
to be caused by nematic CDW modulations. This is con¬ 
firmed by looking at lower dopings (see below). 

There is a second type of nematicity in YBCO, which is 
not associated with CDW modulations. The signature of 
this nematicity is an in-plane anisotropy in the resistivity 
p{T) that is distinct from the anisotropy due to the CuO 
chains. It shows up as an upturn in the anisotropy ratio 
Pa / Pb at low temperature, first detected at low doping 
(p < 0.1) [19]. In particular, it is observed at dopings 
where there are no CDW modulations (Fig. 1) and no 
Fermi-surface reconstruction [13], below p = 0.08. In 
that region, the Nernst anisotropy, vi, / Va, was found to 
track the conductivity anisotropy, Ub / <Ja (= Pa / Pb) [22] 
(see Fig. 2(a)). This second kind of nematicity may be 
associated with spin degrees of freedom, since an in-plane 
anisotropy in the spin fluctuation spectrum measured by 
inelastic neutron scattering emerges spontaneously upon 
cooling [23], in tandem with the rise in pa / Pb [19]. 

In this Article, we report measurements of p(T) and 
v{T) in YBCO performed in the intermediate doping 
range 0.08 < p < 0.12, and use the data to disentan¬ 
gle the CDW-related anisotropy from the nematicity at 
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FIG. 4: Anisotropy difference in the Nernst coefficient of YBCO, defined as D{T) = {va — Vb) / T, and plotted as D{T) — 
D{200 K) vs T / T* in panel (a), and as D(T) — D{200 K) vs T /Txrd in panel (b), where T* is given by the dashed line 
in Fig. 6 and Txrd by the blue line in Fig. 1. Data at p = 0.12,0.13 and 0.15 are from ref. 22. Panel (a) shows that the 
anisotropy at low temperature is not controlled by the pseudogap temperature T* across all dopings. In panel (b), the collapse 
of the three curves at p = 0.12, 0.13 and 0.15 onto a common curve shows that at those dopings the upturn in D(T) at low T 
is correlated with the onset of CDW modulations (at Txrd), in agreement with Fig. 3. By contrast, at p = 0.085 and 0.11, 
the anisotropy is not linked to the CDW order. Indeed, those two curves at lower doping are correlated instead with Tnem, the 
onset of anisotropy in the resistivity (Fig. 1). 


low doping. We define the onset temperature for the lat¬ 
ter nematicity, T^eim track it as a function of doping 
(Fig. 1). We find that Tnem hits the CDW dome (Txrd) 
near its peak and that it scales with the pseudogap en¬ 
ergy EpG- 


II. METHODS 

Single crystals of YBa 2 Cu 30 y (YBCO) were prepared 
by flux growth [24]. Their hole concentration (doping) p 
is determined from the superconducting transition tem¬ 
perature Tc [17], defined as the temperature below which 
the zero-field resistance is zero. A high degree of oxygen 
order was achieved for a pair of samples with p = 0.11 
{y = 6.54, ortho-II order) and with p = 0.12 (y = 6.67, 
ortho-VIII order). Each measurement of anisotropy was 
done on a pair of samples oriented with the long trans¬ 
port direction along the a and b axis of the crystalline 
detwinned structure, respectively. The in-plane electrical 
resistivity of YBCO was measured at dopings p = 0.11 
and 0.12, in each case on a pair of a axis and b axis 
samples (Fig. 5(b)). The Nernst effect - the transverse 
voltage (along y) generated by a longitudinal thermal 
current (along x) in the presence of a perpendicular mag¬ 
netic field (along z) - was measured as described else¬ 
where [22]. The Nernst response of YBCO was mea¬ 
sured at dopings p = 0.085 {H = 16 T) and p = 0.11 
(H = 18 T), in each case on a pair of a axis and b axis 
samples (Figs. 4, 5(c) and 5(d)). The thermal current 


was induced along the a (6) axis and the voltage was 
measured along the transverse b (a) axis for a {b) axis 
samples. In all measurements, the magnetic field was 
applied along the c axis, normal to the Cu 02 planes. 
For the temperature range of data shown in Fig. 5, the 
Nernst coefficient (and its anisotropy) is independent of 
magnetic field strength. 


III. RESULTS 

In Fig. 5(a), we reproduce data from Ando et al. [19] on 
YBCO at two dopings, and define the temperature Tnem 
as the onset of this upturn. We find Tnem = 250 ± 20 K 
and 210 ± 20K at p = 0.062 and p = 0.08, respectively. 
In Fig. 1, we plot Tnem values so obtained on the phase 
diagram of YBCO. 

In Fig. 5(b), we show the anisotropy ratio Pa / Pb vs 
T measured in our own YBCO crystals. At p = 0.11, 
we observe a clear upturn, below Tnem = 150 ± 20 K . In 
Fig. 1, we see that this value is in smooth continuation of 
the Tnem hue at lower doping. Close inspection of earlier 
data [19] at y = 6.70 {p = 0.104) reveals a faint anomaly 
at T ~ 160 K, which our oxygen-ordered samples bring 
out more clearly. With Tnem Txrd at p = 0.11, we 
discover that the nematic line hits the CDW dome near 
its peak (Fig. 1). Going up to p = 0.12, we see no sign 
of any upturn in pa / Pb (Fig. 5(b)), in agreement with 
previous work [19,22] (Fig. 2(b)). 

In order to detect Tnem in a second observable, we mea- 
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FIG. 5: (a, b) Temperature dependence of the anisotropy in the in-plane resistivity of YBCO above T^, measured as the ratio 
of Pa to pb, the electrical resistivities perpendicular {pa) and parallel {pb) to the CuO chains of the orthorhombic structure. 
We dehne the onset temperature for nematicity, Tnem (down-pointing arrows), as the start of the upturn as temperature is 
reduced. Tnem is then plotted in the phase diagram of Fig. 1. (a) Zoom on the data of ref. 19, for two dopings as indicated, 
(b) Anisotropy ratio pa / Pb in our samples, at two dopings as indicated. At p = 0.11 (blue), an upturn is clearly seen, below 
Tnem = 150 ± 20 K (down-pointing arrow). At p = 0.12 (red), no upturn is detected. In the p = 0.11 curve, the drop in pa / pb 
below 80K coincides with the drop in the Hall coefficient at p = 0.11 (ref. 13), due to Fermi-surface reconstruction, (c, d) In¬ 
plane anisotropy difference in the Nernst coefficient Va:{T), plotted as {va — Vb) / T (blue dots), for p = 0.085 {H = 16 T; c) and 
p = 0.11 {H = 18T; d). The data shown here are independent of the magnetic field strength, H. The red dashed lines are a 
fit to the linear dependence at high temperature. The red arrows mark the onset of a large upturn as temperature is reduced, 
at temperatures that fall on the Tnem line of Fig. 1 (red squares). 


sured the Nernst coefficient Vx{T) in YBCO at p = 0.085 
and p = 0.11. In Figs. 5(c) and 5(d) we plot the in-plane 
anisotropy difference, as {ua — I'b) / T vs T. (Note that 
our data at p = 0.085 are in good agreement with data at 
p = 0.078 from ref. 22.) With decreasing temperature, we 
can readily identify the onset of an upturn at 200 ± 20 K 
and 140 ± 20K for p = 0.085 and p = 0.11, respectively. 
These onset temperatures are in excellent agreement with 
the Tnem Hne (Fig. 1), thereby confirming the location of 
this line with a thermo-electric measurement. 

In Fig. 4, we compare the Nernst anisotropy at five 
dopings, from p = 0.085 to p = 0.15. We see that the 
curves split into two groups: those associated with CDW 
modulations, at p = 0.12 and above, whose upturn is 
correlated with Txrd , and those associated with the sec¬ 
ond kind of nematicity, at p = 0.11 and below, which 


do not collapse onto the common curve (whether scaled 
by Txrd or T*), but start their upturn at Tnem- This 
shows that the upturn in (— r'h ) / T is not correlated 
with Txrd, or with T*, across all dopings. We see that 
there are in fact two types of nematicity in YBCO, with 
distinct characteristic temperatures in the phase diagram 
(Fig. 1): one associated with Tnem - seen in both v and 
p - and another associated with Txrd “ seen only in v. 
In the remainder, we focus on the former nematicity, at 

p < 0.11. 

It is instructive to overlay our Tnem data on the cur¬ 
vature map produced by Ando et al. [25], where they 
plot the second derivative of Pq(T) throughout the T — p 
plane. This is done in Fig. 6. Upon cooling from room 
temperature, Pa(T) in YBCO deviates downward from 
its linear-T dependence at high T below a temperature 
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FIG. 6: Map of the curvature in the resistivity of YBCO, 
showing the second temperature derivative of pa{T) in the 
T — y plane, where y is the oxygen content (from ref. 25). 
Red is for positive curvature, blue for negative curvature, and 
white for no curvature. The superconducting transition tem¬ 
perature Tc is shown as solid green circles, with the corre¬ 
sponding p values given on the top. The onset of nematicity 
at Tnem is showu as red symbols (from Fig. 1). Tnem is seen 
to lie right at the inflexion point in the resistivity, Tx, where 
pa{T) goes from downward (blue) to upward (red) curvature. 
The lower boundary of the white region of linear-T resistivity 
at high temperature (upper right corner) is traditionally used 
to define the pseudogap temperature T* . The white dashed 
line is the T* line obtained from Nernst data [22]. (The blue 
band above Tc, of width 20-30 K, is due to paraconductivity, 
i.e. superconducting fluctuations.) 


labeled T* — the upper boundary of the blue region of 
negative curvature. Traditionally called the “pseudogap 
temperature”, T* is also detected as a drop in j T 
(dashed line; from ref. 22). It is roughly the temperature 
below which the NMR Knight shift Ks(T) starts to drop 
in YBCO [26], for example. 

As Pa{T) decreases below T*, it goes through a well- 
defined inflexion point, at Tx — the white line in Fig. 6 
between the blue region of negative curvature at interme¬ 
diate temperatures and the red region of positive curva¬ 
ture at low temperature. We immediately notice that the 
Tnem data points coincide perfectly with the inflexion line 
in Pa(T), i.e. Tnem = Tx, at all dopings. This comes from 
the fact that the upturn in the anisotropy ratio Pa / Pb is 
due to the upward curvature in Pa{T) itself, which starts 
at Tx. This reinforces the fact that there exists a real 
line in the phase diagram that is distinct from the CDW 
dome and separate from T*. 

Note that nematicity was also observed in 
La 2 -a;Sr 3 ;Cu 04 (LSCO) at low doping, with com¬ 
parable values of Tnem [19] I and the region of upward 
curvature in the curvature map of LSCO has a boundary 


similar to that of YBCO [25]. Moreover, the in-plane 
resistivity of the cuprate HgBa 2 Cu 04+5 (Hg-1201) has 
an inflexion point very similar to that found in YBCO 
at Tx [27]. It has been emphasized that PaiT) has a 
Fermi-liquid-like T^ dependence below a temperature 
T** ~ Tx in both YBCO and Hg-1201 [28]. From these 
similarities, we infer that the onset of nematicity at 
Tnem = Tx is likely to be a generic property of hole-doped 
cuprates. 


IV. DISCUSSION 

To summarize, two distinct features of p(T) — inflexion 
point and nematicity — are intimately linked, and remain 
linked at all dopings, along a line that decreases linearly 
withp, the Tnem = Tx line (Figs. 1 and 6). This Tnem line 
has two important properties. The first is that it scales 
with the pseudogap energy. Indeed, as shown in Fig. 7, 
Tnem TpQ, where Tpc is the energy gap extracted from 
an analysis of several normal-state properties of YBCO 
(and other cuprates) above Tc, including Ks{T), Pa{T) 
and the specific heat [29] . This is consistent with analyses 
showing that KgiT) and Pa{T) can be scaled on top of 
each other at all dopings, when plotted vs T / Tx [30,31]. 
The simultaneous rise of Tnem and EpQ with decreasing p 
in YBCO (Fig. 7) is reminiscent of the simultaneous onset 
of nematicity and pseudogap in Bi 2 Sr 2 CaCu 2 08-i-5 (Bi- 
2212) detected in low-temperature STM measurements 
when doping is reduced below p ~ 0.19 [32]. 

The second important property of the Tnem line is that 
it hits the Txrd dome at its peak (Fig. 1). This suggests a 
link between nematicity onset and CDW dome. It points 
to a nematic phase that competes with CDW order, pro¬ 
viding a potential explanation for why CDW order and 
modulations peak at p = 0.12. 

Let us consider two possible scenarios for the nematic¬ 
ity that develops in YBCO below Tnem- The first scenario 
is a nematic phase that breaks the rotational symmetry 
of the Cu02 planes [33,34]. Calculations of transport 
anisotropies applied to YBCO showed that the sign and 
magnitude of the measured db / <Ja and Vb / Va are natural 
consequences of such electron nematic order [35]. Here 
Tnem would Correspond to the onset of anisotropic spin 
correlations, observed with neutrons [23,36]. At low tem¬ 
perature, these correlations turn into short-range spin- 
density-wave (SDW) order [36]. The SDW order is known 
to compete with both superconductivity and CDW or¬ 
der [37], hence it may well be responsible for the fall 
of both Tc and Txrd below p ~ 0.12 (Fig. 1). This 
is analogous to a CDW scenario in which Txrd marks 
the onset of “charge nematicity” [20], and the onset of 
long-range CDW order occurs only at the lower temper¬ 
ature Tnmr [1,10]- Recent calculations of the in-plane 
anisotropy of resistivity show how SDW fluctuations can 
give rise to an upturn in Pa / Pb, while CDW modulations 
would produce a downturn [38]. This could explain the 
two distinct behaviors seen in the resistivity anisotropy 
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FIG. 7: Doping dependence of the energy gap Epq ex¬ 
tracted from five normal-state properties above Tc, as indi¬ 
cated (adapted from ref. 29). All data were obtained on 
YBCO except for ARPES data, measured on Bi-2212 [29]. 
EpG is the pseudogap energy, expressed in K. We compare 
EpG with our temperature Tnem (= 7x) (red symbols, from 
Fig. 1), multiplied by a constant factor 2.4. We observe that 
Tnem ~ EpQ, strong evidence that Tnem is a characteristic 
temperature of the pseudogap phase. 

of YBCO for p < 0.11 and p > 0.12. 

In the second scenario, Tnem rnarks the onset of an en¬ 
hanced nematic susceptibility, with no broken rotational 
symmetry. In solutions of the Hubbard model, it has 
been shown that in the presence of a small orthorhombic 
distortion, of the kind present in YBCO, a large in-plane 
anisotropy of the resistivity appears in the pseudogap 
phase [39,40]. This could explain the upturns seen in 
the transport anisotropy of YBCO at low p and low T, 
and the correlation between Tnem and Epq. In these 
calculations, the key organizing principle of the pseudo¬ 
gap phase is the Widom line, a line in the T — p phase 
diagram where electronic properties of the material are 
predicted to have a distinctive anomaly [41]. For exam¬ 
ple, the charge compressibility peaks at the Widom line 
and both the density of states {e.g. the pseudogap mea¬ 
sured in specific heat) and the spin susceptibility {i.e. the 
Knight shift measured in NMR) go through an inflexion 
point [41]. Given that the onset of nematicity in YBCO 
at Tnem coincides with the inflexion point in Pa(T) (at 
Tx), it is tempting to identify the empirical Tnem = Tx 
line of Fig. 1 with the theoretical Widom line of ref. 41. 

V. SUMMARY 

Nematicity in YBCO is detected as the spontaneous 
appearance of an in-plane anisotropy in the transport 


properties upon cooling, in addition to the background 
anisotropy due to the CuO chains that run along the 
b axis of the orthorhombic structure. Two types of ne¬ 
maticity exist in two different regions of the temperature¬ 
doping phase diagram (Fig. 1). At low doping {p < 0.11), 
the first type is detected in both the resistivity and the 
Nernst coefficient, as parallel upturns in the anisotropy 
ratios, Pa/Pb and v^/va (Fig. 2). At high doping {p > 
0.12), the second type is detected only in the Nernst co¬ 
efficient, and not in the resistivity (Fig. 2). We attribute 
this second nematicity to the short-range CDW modu¬ 
lations detected in YBCO by X-ray diffraction, which 
grow gradually upon cooling below Txrd (Fig. 1), since 
the anisotropy ratio t'b/t'a at p > 0.12 grows in tandem 
with the CDW intensity (Fig. 3). 

We discuss two possible scenarios for the first type of 
nematicity (Tnem line in Fig. 1). In the first scenario, 
rotational symmetry is broken below Tnem as a precursor 
to the onset of SDW modulations at lower temperature, 
which then also break translational symmetry. In the 
second scenario, there is no broken rotational symmetry 
but only an enhanced nematic susceptibility appearing 
gradually below Tnem = Tx (Fig. 6), due to the emergence 
of the pseudogap phase. This could explain why Tnem 
scales with the pseudogap energy Epq (Fig. 7). 
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